The QGP shear viscosity — elusive goal or just 
around the corner? 
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Abstract. With the new viscous hydrodynamic + hadron cascade hybrid code 
VISHNU, a rather precise (0(25%)) extraction of the QGP shear viscosity (?7/s)qgp from 
heavy- ion elliptic flow data is possible if the initial eccentricity of the collision fireball 
is known with < 5% accuracy. At this point, eccentricities from initial state models 
differ by up to 20%, leading to an O(100%) uncertainty for (?7/s)qgp- It is shown 
that a simultaneous comparison of elliptic and triangular flow, t>2 and W3, puts strong 
<^> constraints on initial state models and can largely eliminate the present uncertainty in 

■ (t?/s)qgp- The variation of the differential elliptic flow V2(pr) for identified hadrons 

. between RHIC and LHC energies provides additional tests of the evolution model. 

Prologue — how to measure (?7/s)qgp : Hydrodynamics converts the initial spatial 
deformation of the fireball created in relativistic heavy-ion collisions into final state 
momentum anisotropics. Viscosity degrades the conversion efficiency e x = — >■ 

£ p = \t xx +tw) °f ^ ne nm d; f° r gi ven initial fireball ellipticity e x , the viscous suppression 
of the dynamically generated total momentum anisotropy e p is monotonically related to 
the specific shear viscosity T]/s. The observable most directly related to e p is the total 
charged hadron elliptic flow vf 1 [Tj. Its distribution in p? depends on the chemical 
composition and p^-spectra of the various hadron species; the latter evolve in the 
hadronic stage due to continuously increasing radial flow (and so does V2(pr)), even 
if (as expected at top LHC energy [2]) e p fully saturates in the QGP phase. When (as 
happens at RHIC energies) e p does not reach saturation before hadronization, dissipative 
hadronic dynamics [3] affects not only the distribution of e p over hadron species and pt, 
but even the final value of e p itself, and thus of vf 1 from which we want to extract r]/s. 
To isolate the QGP viscosity (??/s)qgp we therefore need a hybrid code that couples 
viscous hydrodynamics of the QGP to a realistic model of the late hadronic stage, such 
as UrQMD that describes its dynamics microscopically. VISHNU [5 J is such a code. 



Extraction of (t7/s)qgp from 200 A GeV Au+Au collisions at RHIC: The 

left panel in Fig. [TJ shows that such an approach yields a universal dependence of 
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the ellipticity-scaled total charged hadron elliptic flow, v^je x , on the charged hadron 
multiplicity density per overlap area, (l/S)(dN ch /dy), that depends only on (??/s)qgp 
but not on the details of the initial state model that provides e x and S [6]. Pre- 
equilibrium flow and bulk viscous effects on these curves are small [6]. 
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Figure 1. (Color online) Centrality dependence of eccentricity-scaled elliptic flow [6]. 

The QGP viscosity can be extracted from experimental v f 1 data by comparing them 
with these universal curves. The right panels of Fig. [1] show this for MC-Glauber and 
MC-KLN initial state models [6j. In both cases the slope of the data [7j is correctly 
reproduced (not true for ideal nor viscous hydrodynamics with constant r]/s). Due 
to the ~20% larger ellipticity of the MC-KLN fireballs, the magnitude of vf^ exp /e x 
differs between the two models. Consequently, the value of (t]/s)qgv extracted from 
this comparison changes by more than a factor 2 between them. Relative to the initial 
fireball ellipticity all other model uncertainties are negligible. Without constraining e x 
more precisely, (??/s)qgp cannot be determined to better than a factor 2 from elliptic 
flow data alone, irrespective of any other model improvements. Taking the MC-Glauber 
and MC-KLN models to represent a reasonable range of initial ellipticities, Fig. [1] gives 
1 < 47t(?7/s)qgp < 2.5 for temperatures T C <T < 2T C probed at RHIC. 




Figure 2. (Color online) Eccentricity-scaled elliptic flow as function of impact 
parameter for pions, kaons and protons from single-shot and event-by-event ideal fluid 
evolution of fluctuating initial conditions from the MC-Glauber (left) and MC-KLN 
(right) models. 

VISHNU with (v/ s )qgp = ^ for MC-Glauber and ^ for MC-KLN provides an 
excellent description of all aspects of soft (py<1.5GeV) hadron production (p T - 



3 



spectra and differential ^(pr) for all charged hadrons together as well as for individual 
identified species) in 200 A GeV Au+Au collisions at all but the most peripheral collision 
centralities [8]. Such a level of theoretical control is unprecedented. 

Event-by-event hydrodynamics of fluctuating fireballs: In Fig. [1] we evolved 
a smooth averaged initial profile ("single-shot hydrodynamics"). This overestimates 
the conversion efficiency v-ije [HI [10] ■ Fig. [2] shows that event-by-event ideal fluid 
dynamical evolution of fluctuating fireballs reduces v-ije by a few percent [10]. The 
effect is only ~ 5% for pions but larger for heavier hadrons. We expect it to be less in 
viscous hydrodynamics which dynamically dampens large initial fluctuations. A reduced 
conversion efficiency v^/e from event-by-event evolution will reduce the value of (7//s)qgp 
extracted from vf 1 ] based on what we see in ideal fluid dynamics, the downward shift 
for (?]/s)qgp will at most be of order 0.02-0.03. 

Predictions for spectra and flow at the LHC: The successful comprehensive fit 
of spectra and elliptic flow at RHIC [8] allows for tightly constrained LHC predictions. 
Fig. |3] shows such predictions for both pure viscous hydrodynamics VISH2+1 [IT] and 
VISHNU [12J. A straightforward extrapolation with fixed (?7/s)qgp overpredicts the LHC 




Figure 3. (Color online) Total charged hadron elliptic flow as function of centrality 
(VISHNU, left [12 ) and differential elliptic flow for identified hadrons for 20-30% 
centrality (VISH2+1, right [11]) for 200 A GeV Au+Au collisions at RHIC and 
2.76 A TeV Pb+Pb collisions at the LHC. Experimental data are from [13]. 

vf 1 values by 10-15%; a slight increase of (t)/s)qgp from 0.16 to 0.20 (for MC-KLN) 
gives better agreement with the ALICE data [13] . However, at LHC energies v-i becomes 
sensitive to details of the initial shear stress profile [11], and no firm conclusion can be 
drawn yet whether the QGP turns more viscous (i.e. less strongly coupled) at higher 
temperatures. The right panel shows that, at fixed pr < 1 GeV, V2(pr) increases from 
RHIC to LHC for pions but decreases for all heavier hadrons. The similarity at RHIC 
and LHC of u| h (pr) for the sum of all charged hadrons thus appears accidental. 

Constraining initial state models by simultaneous measurement of v 2 and i> 3 : 

While the ellipticities e 2 differ by about 20% between MC-KLN and MC-Glauber models, 
their triangularities £3 (which are entirely due to event-by-event fluctuations) are almost 
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Figure 4. (Color online) pT-diffcrcntial elliptic and triangular flow from viscous 
hydrodynamics for initial eccentricities from the MC-KLN and MC-Glauber models. 

identical [10]. This suggests to use triangular flow V3 (which is almost entirely [10] driven 
by e 3 ) to obtain a model-independent measurement of (?7/s)qgp- Fig- IH shows v^Pt) 
and v^(pt) for deformed Gaussian fireballs with average eccentricities £2 and £3 (with 
random relative angle) taken from the fluctuating Glauber ( "MC-Glauber-like" ) and 
KLN ( "MC-KLN-like" ) models. It demonstrates that a given set of flow data requires 
shear viscosities that differ by a factor 2 to reproduce ^(pr) and but the same shear 
viscosities in both models to reproduce t>3 (pt)- A good fit by both models to ^(pt) 
produces dramatically different curves for v$(pt), and vice versa. The figure illustrates 
the strong discriminating power for such simultaneous studies and gives hope for a much 
more precise extraction of (?7/s)qgp in the near future. 
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